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Abstract 

In recent work jlj we have found identical behaviour for various spin models on "thin" random 
graphs - Feynman diagrams - and the corresponding Bethe lattices. In this paper we observe that the 
ratios of the saddle point equations in the random graph approach are identical to the fixed point (s) 
of the recursion relations which are used to solve the models on the Bethe lattice. The loops in 
the random graphs thus have no influence on the thermodynamic limit for such ferromagnetic spin 
models. 

We consider the correspondence explicitly for Ising and q state Potts models and also note that 
multi-spin interaction models on cacti admit a similar correspondence with a randomised version of 
the cacti graphs which contain loops. 



1 Introduction 

In a series of papers we have found identical behaviour for various spin models on "thin" random 
graphs - Feynman diagrams -and the corresponding Bethe lattices. The key observation for the approach 
of was made in g], where it was noted that spin models on random graphs could be treated as the 
N — > 1 limit of the appropriate hermitean matrix model. The N — > oo limit in such models is already 
familiar from the theory of 2D gravity, where it generates planar random graphs and the expansion in 
1/N is an expansion in the topology of the graphs. The N — ► 1 limit weights all topologies equally so the 
graph/surface correspondence of the N — > oo limit, which is important for 2D gravity and string theory 
applications, is lost. Nonetheless, the generic random graphs of Q (which we denoted as "thin" graphs to 
distinguish them from the "fat" graphs of the planar limit) are still of interest in a statistical mechanical 
context because their locally tree-like structure Q means that one obtains mean field transitions for 
spin models which live on them. They offer one great advantage over genuine tree-like structures such 
as the Bethe || lattice in that they have no boundary, so one is not forced to discard the majority of 
vertices in a simulation or be subject to the gymnastics of considering only points "deep within" the 
lattice analytically. 

One remarkable fact that has emerged from all of the simulations and calculations in [[ij is that the 
thermodynamic behaviour of spin models with ferromagnetic couplings on the Feynman diagrams was 
identical to that on the corresponding Bethe lattice, even down to non-universal features such as the 
critical temperature. The loops did play a role when one considered antiferromagnetic couplings, which 
offered the possibility of frustration and the presence of spin glass rather than antiferromagnetic ordering 
at low temperatures. The similarity was all the more interesting because the methods of solution for 
the two classes of model are, at first sight, completely different. In the random graph models one looks 
at a saddle point equation for the "fields" of the model and the saddle point action, which is exact 
in the n = ^vertices — > oo limit, gives the free energy. All the thermodynamic quantities such as the 
magnetisation, energy, specific heat and susceptibilities may then be calculated from this. For the Bethe 
lattice models on the other hand, the hierarchical structure of the lattice allows one to construct a 
recursion relation for partition functions on sub-trees that is iterated shell by shell to get the behaviour 
at the central point (i.e. "deep within" the lattice). The asymptotic behaviour of this iteration (fixed 
points, limit cycle, chaotic . . . ) then indicates the phase structure of the model || fij and provides a 
fruitful tie-in with the theory of dynamical systems, particularly for multi-spin interaction and frustrated 
models. 

In this paper we show that the previous heuristically observed identity between the Bethe lattice 
and its random graph equivalent is explained by the fact that the ratio of the saddle point equations in 
the random graph model gives exactly the fixed point (s) of the recursion relation in the Bethe lattice 
models. The statement that we have made, and repeated, in various of jlj that the loops are irrelevant 
for ferromagnetic phase transitions on random graphs is thus put on a firmer footing. In essence, the 
content of the saddle point equations for the random graph models is identical to the recursion relations 
of the Bethe lattice, which of course contains no loops by definition. We demonstrate this explicitly for 
the Ising and Potts models and, in order to show that the result is not confined to nearest neighbour 
interactions and the Bethe lattice proper, we also consider a three-site interaction model J|, [|, ||, |§] on 
a cactus graph and its randomised equivalent and find exactly the same equivalence. For simplicity, we 
shall consider 3-regular (<^ 3 ) random graphs and cacti built from triangles throughout this paper, which 
correspond to a Bethe lattice with 3 neighbours for each site, but the results we discuss hold for any 
number of neighbours. 



2 The Ising Model 



The partition function for the Ising model with Hamiltonian 
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(<Jij — ±1) on 3-regular random graphs with 2n vertices may be written as |) 



where iV n is the number of undecorated (no spin) </> 3 graphs with 2n vertices 
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the action is 
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where a, 6 = =b and the propagator is 



where g = exp(2/3) relates the coupling in the action to the inverse temperature /?. 
The saddle point equations dA/dcf>± = may be written, after scaling out A, as 

<t>- - Vg^ + ^l (6) 
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and one finds the solutions 
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l + g ± v /(g+l)(g-3). 
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The upper high temperature paramagnetic solution bifurcates to the low temperature ferromagnetic 
solution in a mean-field transition at g — 3, as can be seen by looking at the magnetisation order 
parameter, which in this notation is 

4>l - 3 

At first sight these solutions, and the method by which they were obtained bear little relation to the 
Bethe lattice calculation Q where the partition function may be expressed in a recursive manner as 

Z = ^2 exp(/io- ) [5„Oo)] 3 (9) 
where the central spin is <7o, the external field is h, we have n "shells" in the Bethe lattice and 

9n(<T0) =^Q„(<7 |S) (10) 

s 

with 

2 

Qn(tr \s) = exp(/3a oSl + h Sl ) TJ Q„_i(si|*«) (11) 



1 In principle the thin graph partition functions are for an annealed ensemble of all graphs with 2n vertices. In practice 
self-averaging appears to ensure that one can consider a single random graph, at least for ferromagnetic spin models such 
as those discussed here. 



2 



reflecting the decomposition of the tree into trunk and branches, s denotes all the spins apart from <jq 
on the jth sub-tree and v^' denotes the spins apart from s\ on the jth branch of the sub-tree. Letting 

9n{+) 

and setting the external field h to zero, we find finally 

*•» ~ ^TP? (13) 

This recursion relation gives a single fixed point in the paramagnetic phase (g < 3) and a pair of stable 
fixed points with one unstable fixed point in the ferromagnetic phase (g > 3). 

Although the phase transition points are identical they have apparently appeared by very different 
paths on the Bethe lattice and on the random graphs. However, a much closer parallel appears between 
equ.@ and equ. (|l3|) when one divides the saddle point equation for </>_ by that for <fr + to get 

0_ _ V9t 2 - + 



which is identical to equ.(|l3]) when we set x n = (f)—/<f>+ and go to the fixed point(s). The loops in 
the random graphs have thus made no contribution to the saddle point equations, which are completely 
equivalent to the fixed point solutions of equ.Jl^). 

Although we have only detailed here the results when the external field is zero, the divided saddle 
point equations are still identical to the Bethe lattice recursion relation at the fixed point in non-zero 
field when the identification x — <f>-/4>+ is made. As a consequence, quantities such as the magnetisation 
are also given by identical formulae on thin graphs and on the Bethe lattice (for a site "deep within" the 
lattice). The Ising model is the simplest decoration that we could consider for our Bethe lattice and <fi 3 
random graphs, but the correspondence between the critical behaviour on the two classes of lattices is 
no fluke and is maintained for other spin models, as we see in the next two sections. 



3 The Potts Model 

The q state Potts model Hamiltonian may be taken to be 

H = P ^(W, -1), (15) 



where the spins <7j now take on q different values. The general formula for the partition function is similar 
to the Ising model, but there are now q fields <f> and the action is given by 



^X>?-c£^-^X>?, (16) 
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where the coupling is related to the inverse temperature by 

1 



(exp(2/?) +(7 -2) 
One finds the low temperature saddle point solutions 



1 - (g - 3)c - y/1 - 2(q - l)c + (g - 5)(g - l)c 2 
h... q -i - 2 

, _ 1 + (g - l)c + 0^= 2(q - l)c +(q-5)(q- l)c 2 



(17) 



(18) 
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(and another low temperature branch with the signs in front of the square roots reversed) which can 
be obtained self-consistently by imposing the observed symmetry breaking pattern on the full action in 
equ.(16) to get 

A = \{q - 1) [1 - c(q - 2)] 4? ~ \{q ~ 1)4? + \tf - \ft ~ <q - 1)# (19) 

with 4> = 4>i...q-i- 4 — 4>q an d A scaled out. Similarly, at high temperature the saddle point solution is 
4>i... q = 4o — 1 ~~ {q ~ l) c - which can be obtained from the effective action 

A = |(l-c( g -l))^-|^. (20) 

Although we started with q fields the final saddle point solution only requires two fields in its solution 
as a consequence of the symmetry breaking pattern in the low temperature phase, just as for the Ising 
model. 

The ratio of the saddle point equations following from equ. ( p"9| ) may be written as 

(21) 



4 c(q - 1)4 2 + (1 - c(q - 2))4 2 

and it is not immediately apparent that this is equal to the recursion relations for the Potts model on a 
Bethe lattice in [|l Of] , which may be written in the form 



1 + x n /c 



(22) 



Xn+1 [exp(2p) + (q-l)x rl 
However, one can rearrange equ.(|2l]) into 

1 + Z2/C (23) 



(q - l)z 2 + exp(2^) 

where z — 4/4>- If we now make the identification z 1 = x we recover the fixed point of the recursion 
relation in equ.(]2^). We have thus found, just as for the Ising model, (which is, after all, a q = 2 state 
Potts model) that the fixed points of the Bethe lattice recursion relations and the saddle point equations 
on thin graphs are identical. 



4 Multi-site interaction models and Loopy Cacti 

A Husimi Tree (or "cactus" graph) is a tree-like structure in the large rather like the Bethe lattice, but 
the individual components are polygons joined at their vertices, rather than simply vertices and links. It 
is characterised by the number of polygons articulated at each vertex, 7. An cactus with 7 = 2 is shown 
in Fig. 1 along with the Bethe lattice for which it forms a sort of medial graph. For 7 > 2 one does not 
have such a direct correspondence but the hierarchical structure still, of course, remains which allows 
one to deploy similar techniques to the Bethe lattice solution, using recursion relations that operate shell 
by shell. Since the individual units are now more complicated than vertices, it is possible to put models 
with multi-site interactions on such trees and investigate their phase behaviour |7| Writing down 

a thin graph style action for a such models corresponds to allowing ( predominantly large) loops in the 
branches of the cacti, in the same manner as the Feynman diagram expansion introduces large loops into 
a locally tree like structure. Given the results of the previous sections relating the critical behaviour on 
the Bethe lattice and thin graphs one might expect that the loopy cacti would display exactly the same 
behaviour as the loop-less cacti. 

As an example we take the following multi-spin interaction model (a special case of the model solved 
by Wu and Wu || on the Kagome lattice) 

H = j3{J^cna j a k + h^u i ) (24) 

A 1 
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where the three spin sum is over triangles A. If one considers the possible spin configurations on the 
triangle of the Husimi tree, namely (+ + +), (— — — ), (+ + — ), (+ — —), the appropriate action can 
be written for the lattice which looks locally like Fig.l (but has large loops of triangles) by thinking of 
the spins as residing on the links of the underlying 3 graph. One then has "ghost" links joining "ghost" 
vertices at the centre of each of the triangles, which gives 

A = f 4 + \P- ~ - - fV + - (25) 

where /i = exp(2/3/i), z = exp(2/3J 3 ) and we have again pre-emptively scaled out the vertex counting factor 
A for simplicity. Note that an edge on which a positive spin resides, represented by the "propagator" for 
4> + , picks up a factor of 1/fi to avoid double counting the external field contribution where two triangles 
join. If we now consider the saddle point equations for the model 

4>+ = z[i 2 (j) 2 + + 2[i(j) + (f)- + z<j)_ 

<t>- = n 2 cf> 2 + + 2zfi<p + <p- + (jP_ (26) 

we find that their ratio is again exactly the recursion relation derived for the model on a loop-less cactus 
(Husimi tree) in ||. 

The Hamiltonian may be extended to include nearest neighbour interactions as well 

H = 0(J 3 cndjCTk + J2^2 + h'^cj l ). (27) 

A ij i 

and similar reasoning to the pure 3-spin interaction case above leads to the action 

A = v (f # + \^ - M 2 ^- - z^U+ - ^0 3 _) , (28) 

for the model on loopy cacti, where the additional parameter v — exp(2/3J2). The overall factor of v is 
irrelevant for solving the saddle point equations and we find that their ratio is identical to the recursion 
relation presented for the model on a loop-less cactus in || 



zn 2 v 2 4> 2 + + 2^4> + <j)- + z(f> 2 _ 
{i 2 <p 2 + + 2zn(j> + 4>- + v 2 (t) 2 _ 



(29) 



The loops in the cacti have thus not altered the critical behaviour of the models from the loop-less case, 
at least for ferromagnetic couplings, where one is looking for fixed points of the recursion relations rather 
than more exotic behaviour. 



5 Discussion 

We have seen that for all the models considered - Ising, Potts, multi-site interaction, and mixed inter- 
action, the randomisation of the Bethe lattice or Husimi tree by the introduction of loops has no effect 
on the thermodynamic limit for ferromagnetic phase transitions. This irrelevance of loops is also to be 
expected on graph theoretic grounds || , for it is known that the loop distribution on a regular random 
graph is very strongly weighted to large loops This provides a justification for calling the thin graphs 
locally tree-like. 

The dynamical systems aspects of the recursive approach on Bethe lattices and Husimi trees are lost 
in the ratios of the saddle point equations when one moves to randomised structures, so phases with 
staggered order (in particular, antiferromagnetic order) which appear as cycles in the solution of the 
recursion relations are also absent. This is intuitively reasonable, because the (f> 3 graphs contain both 

2 More precisely, on 3 regular graphs the number of i-cycles, that is loops with i edges, is asymptotically a Poisson random 
variable with mean 2' /(2i). 
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odd and even loops, so the shell by shell alternation of spin values that characterises antiferromagnetic 
order on the Bethe lattice proper would lead to frustration via the loops. 

Indeed, there is strong evidence 0, ^ both analytical and numerical that the Ising model with purely 
antiferromagnetic couplings displays a spin glass phase at sufficiently low temperature on 4> 3 graphs rather 
than the antiferromagnetic phase of the Bethe lattice. It would be extremely interesting to investigate 
some of the more esoteric dynamical phenomena in the recursion relations of the multi-site interaction 
models with antiferromagnetic couplings that are highlighted in ||, [| in order to see what, if any, 
phenomena they corresponded to in the randomised/loopy case. 
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Figure 1: A section of a Husimi tree (clotted) composed of triangles with 7 = 2. The underlying Bethe 
lattice (solid) is also shown, along with the "ghost" nodes at the centre of each triangle. 
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